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Q\ • ABSTRACT 



We have investigated the structural and dynamical properties of triaxial stellar 
systems whose surface brightness profiles follow the r^/" luminosity law - extending 
\ the analysis of Ciotti (1991) who explored the properties of spherical r-"^/" systems. 

A new analytical expression that accurately reproduces the spatial (i.e. deprojected) 
luminosity density profiles (error less than 0.1%) is presented for detailed modelling 
of the Sersic family of luminosity profiles. We evaluate both the symmetric and the 
' non-axisymmetric components of the gravitational potential and force and compute 

I the torques as a function of position. For a given triaxiality, stellar systems with 

. smaller values of n have a greater non-axisymmetric gravitational field component. 

We also explore the strength of the non-axisymmetric forces produced by bulges with 
differing n and triaxiality on systems having a range of bulge-to-disc ratios. The 
increasing disc-to-bulge ratio with increasing galaxy type (decreasing n) is found to 
O ' heavily reduce the amplitude of the non-axisymmetric terms, and therefore reduce 

the possibility that triaxial bulges in late-type systems may be the mechanism or 
^2 ' perturbation for non-symmetric structures in the disc. 

Using seeing-convolved r^/"-bulge plus exponential-disc fits to the K-band data 
from a sample of 80 nearby disc galaxies, we probe the relations between galaxy type, 
Sersic index n and the bulge-to-disc luminosity ratio. These relations are shown to 
. , be primarily a consequence of the relation between n and the total bulge luminosity. 

. In the K-band, the trend of decreasing bulge-to-disc luminosity ratio along the spiral 

Hubble sequence is predominantly, although not entirely, a consequence of the change 
in the total bulge luminosity; the trend between the total disc luminosity and Hubble 
type is much weaker. 

Key words: celestial mechanics: stellar dynamics - galaxies: structure of - galax- 
ies: photometry - galaxies: elliptical and lenticular, cD- galaxies: spiral - galaxies: 
kinematics and dynamics 



in 
oo 



o 



1 INTRODUCTION 

As the quality of photometric data has improved over the 
years (largely due to the use of CCDs), the applicabil- 
ity of a fitting-function which can account for variations 
in the curvature of a light profile has been demonstrated 
for elliptical galaxies (Capaccioli 1987, 1989; Davies et al. 
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1988; Caon, CapaccioUi & D'Onofrio 1993; Young & Cur- 
rie 1994; Graham et al. 1996), and for the bulges of spi- 
ral galaxies (Andredakis, Peletier & Balcells 1995; Seigar 
& James 1998; Moriondo, Giovanardi & Hunt 1998; Khos- 
roshahi, Wadadekar & Kembhavi 2000; Prieto et al. 2001; 
Graham 2001; MoUenhoff & Heidt 2001). These systems are 
not universally described with either an exponential profile 
or an r^''* law (de Vaucouleurs 1948, 1959), but rather a 
continuous range of light profile shapes exist which are well 
described by the Sersic (1968) r^^" model. 
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In ellipticals, the shape parameter n from the Sersic 
model is strongly correlated with the other global proper- 
ties derived independently of the r^''" model, such as: total 
luminosity and effective radius (Caon et al. 1993; Young & 
Currie 1994, 1995; Jerjen & Bingelh 1997; Trujillo, Graham 
& Caon 2001), central velocity dispersion (Graham, Trujillo 
& Caon 2001) and also central supermassive black hole mass 
(Graham et al. 2001). Additionaly, the spiral Ifubble type 
has been shown to correlate with the bulge index n such 
that early-type Spiral galaxy bulges have larger values of n 
than late-type spiral galaxy bulges (Andredakis et al. 1995; 
Graham 2001). This correlation arises from the fact that the 
index n is well correlated with the bulge-to-disc luminos- 
ity ratio (B/D; see, e.g. Simien & de Vaucouleurs 1986) and 
this is one of the parameters used to establish morphological 
type (Sandage 1961). 

Given the abundance of observational work and papers 
now using the Sersic model, it seems timely that a theoret- 
ical study is performed on realistic, analytical models fol- 
lowing the r^''" law. Structural and dynamical properties of 
isotropic, spherical galaxies following r^^" models have al- 
ready been studied in detail in the insightful paper of Ciotti 
(1991). However, most elliptical galaxies and bulges of spiral 
galaxies are known to be non-spherical objects. Typically, 
the mass models which have been used for the study of triax- 
ial galaxies have followed analytical expressions which were 
selected to reproduce the properties of the de Vaucouleurs 
r^/" profile (e.g. Jaffe 1983; Hernquist 1990; Dehnen 1993), 
or more recently the modified Hubble law (Chakraborty & 
Thakur 2000). For that reason, previous studies based on 
these kinds of analytical models, although certainly useful, 
are however unable to probe the full range of properties 
which are now observed in real galaxies. Consequently, it is 
of importance to know how much room for improvement ex- 
its in the study of triaxial objects following the r^''" family 
of profiles. 

Due to the fact that the observed r^''" luminosity pro- 
files cannot be deprojected to yield analytical expression 
for the spatial density^ the r^''" law has been considered 
less useful for studies of detailed modelling. An analytical 
representation (approximation) for the mass density profiles 
which accurately reproduces the observed r^^" luminosity 
profiles would be of great interest for simulations of real 
galaxies. We have therefore derived just such an analytical 
expression for the mass density profiles of the Sersic family of 
models. Our approximation surpasses the accuracy of both 
the Dehnen models for representing the specific r^''* profile 
and also their extension to the double power-law models of 
Zhao (1997). 

In this paper we present a detailed study of how the 
physical properties of triaxial stellar systems change as a 
function of the index n. An accurate analytical expres- 
sion for modelling the spatial density is presented in Sec- 
tion 2. In Section 3 we explore the axisymmetric and the 
non-axisymmetric components of the potential, forces and 
torques associated with a Sersic light distribution. Finally, 
by using literature available K-band observations of a sam- 



t Recently, Mazure & Capelato (2002) have provided an exact 
solution for this, and other related spatial properties, in terms of 
the Meijer G functions when the Sersic index n is an integer. 



pie of 80 spiral galaxies, the physical basis to the n-T (or 
n-B/D) relation is investigated in Section 4. 



2 THE r^/" MODEL 

The projected, elliptically symmetric Sersic r^^" intensity 
distribution 7(r) can be written in terms of the projected, 
elliptical radial coordinate ^ (see details in Trujillo et al. 
2001) such that: 



/(C) = /(0)e-^"«/'-=''''"\ 



(1) 



where 7(0) is the central intensity, and is the effective 
radius of the projected major-axis. Curves of constant ^ on 
the plane of the sky are the isophotes. The quantity 6„ is a 
function of the shape parameter n, and is chosen so that the 
effective radius encloses half of the total luminosity. The ex- 
act value is derived from r(2n)=27(2n, b„), where r(a) and 
7(0,2;) are the gamma function and the incomplete gamma 
function respectively (Abramowitz & Stegun 1964, p. 260). 
The index n increases monotonically with the central lumi- 
nosity concentration of the surface brightness distribution 
(Trujillo, Graham, & Caon 2001). 

The total projected luminosity L associated with this 
model is given by 



2-Kn 



r(2n). 



(2) 



where e is the ellipticity of the isophotes. For a homologous 
triaxial ellipsoid, the spatial (deprojected) luminosity den- 
sity profile v{C,) can be obtained by an Abel integral equation 
(Stark 1977): 



rl/2 
TT 



(3) 



where f^^'^ is a constant that depends on the three- 
dimensional spatial orientation of the object (Varela, 
Munoz-Tunoz & Simonneau 1996; Simonneau, Varela & 
Munoz-Tunoz 1998) and C, parametrizes the ellipsoids of 
constant volume brightness, f^^'^ equals 1 when the proper 
axis frame of the object has the same orientation as the ob- 
server axis frame (i.e. when the Euler angles between the 
two frames equal zero). 



2.1 Mass density profiles 

Assume a triaxial object whose mass is stratified over ellip- 
soids with axis ratios a:b:c (a>b>c) and the x- (z-) is the 
long (short) axis (see Fig. 1). The symmetry of the problem 
motivates us to work with ellipsoidal coordinates where: 

X = sin 7/; cos 
y — a( sin tp sin 6 

z = P( cos ip (4) 

and where a=b/a and /3=c/a. The mass models considered 
in this study are the triaxial generalizations of the spheri- 
cal models discussed in detail by Ciotti (1991). The math- 
ematical singularities present in Eq. 3 were considered and 
solved by Simonneau & Prada (1999, Eq. 16). Substituting 
Eq. 1 into Eq. 3, letting ^ = ( cosh tp, and multiplying by the 
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Figure 2. a) The dimensionless mass density profiles (see Section 
2.1) for the values of n=0.5, 1, 2, 4 and 10. b) the relative error 
between the the analytical approximation proposed in Eq. 7 and 
the exact solution are shown for the previous values of n. 



Figure 1. A surface of constant density for the triaxial ellipsoid 
described in Eq. (5) and (6). 



mass-to-light ratio T =M/L, we obtain a similar expression 
to the one found by these authors: 



e ""1,^— j (<^coshv3)'/"~'dv',(5) 



with 



\a 



(6) 



The dimensionless mass density profiles p((") = r'i/Mp{C,), 
where M is the total mass, are shown for different values of 
n in Fig. 2a. It should be noted that the inner density pro- 
file decreases more slowly with increasing radius for systems 
having lower values of n. 

The mass density profiles of the r^''" family (Eq. 5) can 
be extremely well approximated by the analytical expres- 



^-''(^) = v:^ ^ — T^W) — ' ^ ^ 

where h = C,/re, C{h) = /ii(log/i)^ + h2 log/i + hy, and Kv 
is the !/th-order modified Bessel function of the third kind 
(Abramowitz & Stegun 1964, p. 374). In the Appendix A we 
show the values of the parameters (j/,p,hi,h2,h3) as function 
of the index n. This approximation contains two exact cases: 
n=0.5 and n=l, and provides relative error less than 0.1% 
for the rest of the cases (Fig 2b) in the radial range 10"'^ < 
C/re < 10"^. This approximation surpasses (by a factor of 
10^-10'') the expression presented in Lima Neto, Gerbal & 
Marquez (1999). 



3 NON-AXISYMMETRIC PERTURBATIONS 



DUE TO A TRIAXIAL 



STRUCTURE 



For three different triaxiality mass distributions: a) spherical 
(q=/3=1); b) moderately triaxial (q:=0.75, /3=0.5); c) highly 
triaxial (a=0.5, /9=0.25), we have explored, in detail, the 



non-axisymmetric gravitational field over the z=0 plane (i.e. 
the disc plane when studying spiral galaxies) . 



3.1 Non-spherical component of the gravitational 
potential in the plane z=0 

We evaluate this quantity by calculating: 



G(r) 



$2(r) 



(8) 



where ^2{r) and ^o{r) are the m=2 and m=0 component 
of the gravitational potential, such that the nth-order term 
$m(r) is evaluated from the gravitational potential on the 
z=0 plane $(r, 6) by using the Fourier decomposition (see, 
e.g. Combes & Sanders 1981). Gravitational potential and 
gravitational force expressions are shown on Appendix B. 

The profiles of G(r) for different triaxialities and values 
of n are shown in Fig. 3. As it is expected, as the triaxial- 
ity increases the non-spherical nature of the gravitational 
field increases. Also, we highlight that for a given triaxial- 
ity, smaller values of n (i.e less concentrated mass distri- 
bution) give greater non-spherical gravitational fields. The 
maximum non-axisymmetrical behavior of the potential is 
obtained at radial distances less than 2 r^. This radial dis- 
tance is also a function of the index n, decreasing a n in- 
creases and remains quite independent of the triaxiality of 
the object. For a moderately triaxial object with n=l, the 
non-axisymmetrical component of the potential can vary 
some 6% between r=0 and r=2re, and varies some 15% for 
our highly triaxial model. 

For an n=l model, and starting from our moderately 
triaxial case (a=0.75, /9=0.50), we increased the value of 13 
to 0.75. The results are shown in Fig. 3c and reveal that 
G(r) varied only mildly. This figure shows that the non- 
axisymmetric effect (along the radial distance) in the z=0 
plane is mainly due to how the mass of the bulge is dis- 
tributed in this plane. 
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3.2 Non— spherical component of radial 
gravitational forces in the z=0 plane 

The non-spherical component of the radial gravitational 
forces in the z=0 can be estimated by: 



N{r) 



d<^2{r)/dr 



(9) 



In Fig. 3 the N(r) profiles (Eq. 9) are evaluated for the 
same cases as was the G(r) profiles. A remarkable point is 
that N(r) reaches its maximum value in the radial range 
2 re<r<4 r^. For a spiral bulge structure, this means that 
the most important non-axisymmetric effects take place in 
a zone which is dominated by the disc. As with the G(r) 
parameter, stronger distortions occur as the triaxiality in- 
creases and the index n decreases. The mechanism which 
controls this distortion is basically determined by the mass 
distribution in the z=0 plane (Fig 3c). 

It is noted that the relative (i.e. % change) non- 
axisymmetric effects on the radial forces are larger than 
the relative distortion on the potential. As an example, 
for a moderately triaxial structure with n=l the non- 
axisymmetric component of the radial forces can reach 8%. 












I \ 





3.3 Torques on z=0 plane 

The torques provoked by the triaxial structures along the 
angular coordinate are evaluated around the circle of radius 
Tmax where the maximum non-axisymmetrical distortion of 
the radial forces is produced (i.e. at the peak of the N(r) 
profile). Given the gravitational potential <l>(r, 6) in the z—0 
plane, we have at the radius Vmax'- 



Ft {r max) 
Fpi {rmax ) 



(10) 



where FTifmax) = [d<^{rmax,6) /d6]/rmax represents the 
amplitude of the tangential force along the angular coor- 
dinate at radius rmax, and FR{rmax) = {d^{rmax,0)/dr) is 
the radial force at this radius. Due to the symmetry of the 
ellipsoid, the values of Il{6) need only be plotted for one 
quadrant in the z=0 plane; we use 0°< 6' <90° (Fig. 3). De- 
pending on the quadrant, Il{6) is either negative or positive 
because the sign of the tangential force changes from quad- 
rant to quadrant. The maximum torque around a circle of 
radius rmax depends on the triaxiality of the object. As the 
triaxiality increases the maximum torque tends to be closer 
to the major axis - as one would expect. The position of 
this peak is quite independent of the value of n. 

The absolute value of the torque for any given triaxi- 
ality increases as n decreases. For our highly triaxial bulge, 
U{e) ranges from 0.17 (n=10) to 0.24 (n=l), which would 
be considered a "bar strength" class of 2 in the classification 
scheme of Buta & Block (2001). In the case of our moder- 
ately triaxial object, the maximum absolute value of Ti{9) 
ranges between 0.06 and 0.09. These values correspond to a 
"bar strength" class of 1. Thus, even a moderately triaxial 
bulge is capable of provoking non-negligible torques on a 
disk - that is to say, a bar is not neccessarily required. A 
detailed study separating the torque contribution from both 
bars and bulges would of course be of interest, and it is our 
intention to add a range of bar potentials to our models in 
the future. 



Figure 3. The parameters G(r), N(r) and Tl{9) are shown for 
different values of n and triaxiality: a) First Row: G(r), N(r) and 
n(9) for a moderately triaxial object and different n; b) Second 
Row: G(r), N(r) and 11(6) for a highly triaxial structure; c) Third 
Row: G(r), N(r) and 11(6) for three moderately triaxial objects 
with n=l and the same axis ratio along the y and x axis. 



As with the previous parameters, for the range of triaxi- 
alities investigated and a given n, varying the mass distribu- 
tion along the z-axis (i.e. varying the triaxiality parameter 
/3) only results in a slight change to Il{9) (see Fig. 3c). For 
a spherical distribution all above parameters are 0. 



4 LINKING THEORY AND OBSERVATIONS: 
THE CONNECTION BETWEEN n AND THE 
B/D LUMINOSITY RATIO 

In the previous section we have seen how the non- 
axisymmetrical effects (in the z=0 plane) from a triaxial 
bulge increase as n decreases. Taken with the correlation 
between n and galaxy type (Andredakis, Peletier & Balcells 
1995) shown in Fig. 4, this invokes the natural question: 
How, if at all, are the structural properties of the bulges 
(i.e. n) related (that is, cause and effect) with the non- 
axisymmetrical components (i.e. arms) observed in the disc? 
The results obtained in the previous sections were evalu- 
ated without any mention to the relative mass of the bulge 
and disc. It turns out that the axisymmetrical mass dis- 
tribution of the disc causes a strong softening of the non- 
axisymmetrical perturbation caused by the non-sphericity 
of the bulge. The degree of "smoothing" is an increasing 
function of the D/B ratio. 

Fig. 5 shows the N(r) profile for a moderately triaxial 
bulge with n = 1 and B/D=0.1 and 0.01, and for a bulge 
with n = 4 and B/D^l.O and 0.1 (Fig. 4). N(r) was evalu- 
ated here assuming the disc follows an exponential surface 
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Figure 4. The best-fitting bulge index n is plotted against the 
B/D luminosity ratio (left) and versus the galaxy morphologi- 
cal type (right). The galaxies come from the samples of de Jong 
(1996) (solid points) and Andredakis et al. (1995) (open dia- 
monds). See text for details. 
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Figure 5. The N(r) profiles for a moderately triaxial bulge with 
n=l and 4, and B/D=l, 0.1, 0.01 are shown. 



brightness distribution. The ratio between the length scale 
of the disc and the effective radius of the bulge is assumed to 
be constant and with a value of h/r^=^ 0. The expressions 
for the potential and the radial force of these structures can 
be found in Binney & Tremaine (1987, p. 77 and 78). 

Fig. 5 illustrates that for B/D luminosity ratios typi- 
cal of real galaxies, the non-axisymmetrical effects on the 
disc largely disappear (< 3%). For B / D=l, the values of 
the N(r) profile remain basically unchanged to that seen in 
Fig. 3, but for B /D=Q.l these values decrease approximately 
by a factor 2, and for i3/D=0.01 this factor is bigger than 
10. Thus, although the non-radial effects on the z=0 plane 
increase as n decreases, the smoothing effects of the increas- 
ingly dominant disc are stronger. Bulges with small values 
of n are unable to produce significant non-axisymmetrical 
effects on a massive disc. 



t Although there is a range of bulge-to-disc size ratios, a median 
value for h/ve in the if -band is 5 (Graham & Prieto 1999). 



4.1 Why does the n—Type (or n—B/D) relation 
exist? 

To explore the connection between bulges and discs in spi- 
ral galaxies (see, e.g. Fuchs 2000), we have used the data 
from two independent samples of galaxies observed in the 
K-band. The K-band provides a good tracer of the mass 
due to the near absence of dust extinction and the reduced 
biasing effect of a few per cent (in mass) of young stars. 
We used the data from Andredakis et al. (1995) and the 
structural analysis of the de Jong (1996) data performed 
by Graham (2001). Both studies were done fitting a seeing 
convolved Sersic law to the spiral galaxy bulges. In both 
samples we have removed those objects which contained a 
clear bar structure, leaving a total of 28 objects from An- 
dredakis et al. (1995) and 52 objects from Graham (2001). 
The relations present in Fig. 4 between n and the B/D lu- 
minosity ratio, and n versus the morphological type T have 
a Spearman rank-order correlation coefficient of Ts = 0.77 
and Ta — —0.73, respectively, for the combined sample. 

Andredakis et al. (1995) suggest "although other possi- 
bilities cannot be excluded, the most straightforward expla- 
nation for this trend is that the presence of the disc affects 
the density distribution of the bulge in such a way as to make 
the bulge profile steeper in the outer parts. One mechanism 
to produce such an effect might be that a stronger disc trun- 
cates the bulge, forcing its profile to become exponential". 
Following this line of thought, via coUisionless N-body sim- 
ulations, Andredakis (1998) studied the adiabatic growth of 
the disc onto an existing r^''* spheroid. He found that the 
disc potential modifies the bulge surface brightness profile, 
lowering the index n. This decrease was larger with more 
massive and more compact discs. This mechanism, however, 
saturated at around n—2 and exponential bulges could not 
be produced. 

We believe that this line of reasoning is not the most 
appropiate explanation for the relation between n and B/D. 
Firstly, we find that the index n is not only well correlated 
with the luminous B/D ratio, but is equally well corre- 
lated {vs — —0.75) with bulge luminosity MK^ibulge) (Fig. 
6a). Additionally, the correlation between n and disc lu- 
minosity MKridisc) is relatively poor (r^ = —0.53). Sec- 
ondly, MKj.{bulge) is more strongly correlated (vs — —0.86) 
with the B/D ratio than MK^idisc) and the B/D ratio 
{ts = —0.50) (Fig. 6b). Hence, it is variations in the bulge 
which are predominantly responsible for variations in the 
B/D ratio. 

These above two correlations seem to indicate that n 
may be related directly with the properties of the bulge 
rather than with the combined B/D ratio. Consequently, 
as n is correlated with the total bulge luminosity, the cor- 
relation between n and B/D is a result of the more fun- 
damental correlation between Mkt {bulge) and B/D. That 
is, it is not the relative increase in disc-to-bulge luminosity 
which produces bulges with smaller values of n, but simply 
that bulges with larger values of n are more luminous (or 
vice-versa) and this produces the correlation between n and 
the B/D luminosity ratio. 

Favouring this argument, we note that among the El- 
liptical galaxies (without the need to invoke any disc) there 
exists a strong correlation (Pearson's r=-0.82; Graham, Tru- 
jillo & Caon, 2001) between n and the total luminosity of 
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Figure 6. a) The absolute K-band magnitude of the bulge (top 
panel) and the disc (bottom panel) are plotted versus the index 
n. The galaxies come from the samples of de Jong (1996) (solid 
points) and Andredakis et al. (1995) (open diamonds). See text 
for details, b) The absolute K— band magnitude of both bulge and 
disc are shown as a function of the B/ D ratio. The galaxies come 
from the samples of de Jong (1996) (solid points) and Andredakis 
et al. (1995) (open diamonds). See text for details. 



these objects. The index n of pressure supported stellar sys- 
tems are related with the total luminosity of these struc- 
tures. In agreement with this, Aguerri, Balcells & Peletier 
(2001) have found (using coUisionless N-body simulations) 
that the bulges of late-type galaxies can increase their n 
values via dense satellite accretions where the new value of 
n is found to be proportional to the devoured satellite mass. 



4.2 (bulge) versus B/D for classifying 

morphological types 

Due to the strong correlation between the B/D luminos- 
ity ratio and Mkt (bulge), it might be of interest to ask 
which one of these quantities is preferred to establish the 
morphological type T of a galaxy]^ Working from B-band 
images (which is good for observing the young star pop- 
ulation, and consequently the spiral arm structure, which 
is one of the basic criteria to the Hubble galaxy classifica- 
tion), Simien & de Vaucouleurs (1986) fitted r^^^ profiles 
and exponential discs to a sample of 64 spiral galaxies and 
34 SO type galaxies. They presented a good correlation be- 
tween the bulge-to-disc luminosity and T, but not between 
Mbt (bulge) and T. Consequently, their B-band observa- 
tions suggested that the B/D luminosity ratio was preferred 
to Mb-t (bulge) for establishing the morphological type T. In 
Fig. 7 we show the relation between the B/D luminosity ra- 
tio and T (vs = —0.65) and between MKT(bulge) and T 
(va = 0.67). Thus, from K-band observations, and fitting 
r^^" bulge profile models, the change in the luminous mass 
of the bulge along the Hubble sequence appears equally as 
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Figure 7. The K-band B/D ratio and the absolute K-band mag- 
nitude of the bulge are shown as a function of the morphological 
type. The galaxies come from the samples of de Jong (1996) (solid 
points) and Andredakis et al. (1995) (open diamonds). See text 
for details. 



important as the combined change in the bulge and disc lu- 
minosityPU It would then follow that the luminous mass of 
the bulge (i.e. AIkj, (bulge)) is related with the spiral arm 
structure. 



5 CONCLUSIONS 

The main results of this work are the following: 

a) We have generalised the analysis of the physical prop- 
erties of spherical stellar systems following the r^^" luminos- 
ity law to a homologous triaxial distribution. The density 
distribution, potential, forces and torques are evaluated and 
compared with the spherical case when applicable (Ciotti 
1991). An extremely accurate analytical approximation (rel- 
ative error less than 0.1%) for the mass density profile is 
provided. 

b) We derive an exact expression showing how the cen- 
tral potential decreases as triaxiality increases. We also show 
that for a fixed triaxiality, as the index n decreases the non- 
axisymmetrical effects in the z=0 plane increase. Even for a 
moderately triaxial object, the non-axisymmetrical compo- 
nent of the potential and the radial forces are not negligible 
for small values of n. These components can range from 6 
to 8%, respectively, compared to the value of the spherical 
component. For our highly triaxial model, they can range 
over some 20%. 

c) The non-axisymmetrical effects in the disc plane due 
to the bulge structure are strongly reduced when an axisym- 
metrical disc mass is added. For this reason, bulges with 
smaller values of n appear unlikely to produce any signifi- 
cant non-axisymmetrical effect on their disc, which is typi- 
cally 10 to 100 more times more massive than the bulge. In 
this regard, the B/D mass ratio and the triaxiality of the 



§ We refer here to the morphological type established on the ba- 
sis of B-band observations. Infrared images have shown that the 
appearence of galaxies can be substantially different (Block et al. 
1999). 



^ We were able to repeat this analysis using B-band data from de 
Jong (1996) (excluding barred galaxies), which gave = —0.64 
between the B/D luminosity ratio and T, but only = 0.54 
between Mbj, (bulge) and T. 
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bulge are more important, that is, can dominate over the 
effects of small n. 

d) The correlation found between n and the B /D lumi- 
nosity ratio found in spiral galaxies is explained here not as a 
consequence of the interplay between the bulge and the disc, 
but due to the strong correlation between n and Mt (bulge), 
and between Mrihulge) and B/D. Also, K-band data do 
not support the idea that the B/D luminosity ratio can be 
preferred over Mt [bulge) as an indicator to establish galaxy 
morphological type (T). Both parameters present equally 
good correlations with galaxy type T. 



REFERENCES 

Abramowitz M., Stegun I., 1964, Handbook of Mathematical 

Functions. Dover, New York 
Aguerri, J.A.L., Balcells M., Peletier R., 2001, A&A, 367, 428 
Andredakis Y.C., Peletier R., Balcells M., 1995, MNRAS, 275, 

874 

Andredakis, Y.C., 1998, MNRAS, 295, 72 

Binney, J., Tremaine, S.: 1987, Galactic Dynamics, Princeton Uni- 
versity Press, Princenton 
Buta R., & Block D.L., 2001, ApJ, 550, 243 

Block, D.L., Puerari, I., Frogel J. A., Eskridge, P.B., Stockton A., 

Burkhard F., 1999, Astropli. and Space Sci, 269-270, 5 
Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 1013 
Capaccioli M., 1987, in Structure and Dynamics of Elliptical 

Galaxies, lAU Symp. 127, Reidel, Dordrecht, p. 47 
Capaccioh M., 1989, in The World of Galaxies, ed. H. G. Corwin, 

L. Bottinelli (Berlin: Springer- Verlag), p. 208 
Chakraborty D. K. & Thakur P., 2000, MNRAS, 318, 1273 
Chandrasekhar, S. 1969, Ellipsoidal Figures of Equilibrium (New 

York: Dover) 
Ciotti, L. 1991, A&A, 249, 99 
Combes, F., & Sanders, R.H. 1981, A&A, 96, 164 
Davies J., Philips S., Cawson M., Disney M., Kibblewhite E., 

1988, MNRAS, 232, 239 
Dehnen, W., 1993, MNRAS, 265, 250 
de Jong, R.S., 1996, A&AS, 118, 557 
de Vaucouleurs G., 1948, Ann. Astrophys, 11, 247 
de Vaucouleurs G., 1959, Hand. Phys., Vol. 53, p. 311 
Fuchs, B., 2000, in: F. Combes and G. Mamon (eds.). Galaxy 
Dynamics: from the Early Universe to the Present, ASP Con- 
ference Series, Vol. 197 
Graham A.W., Lauer T., CoUes M.M., Postman M., 1996, ApJ, 
465, 534 

Graham A.W., Prieto M., 1999, ApJ, 524, L23 
Graham A.W., 2001, AJ, 121, 820 

Graham A.W., Trujillo, I., Caon, N., 2001, AJ, 122, 1707 
Graham A.W., Erwin, P., Caon, N., Trujillo, I., 2001, ApJ, 563, 
Lll 

Hernquist, L., 1990, ApJ, 356, 359 
Jaffe, W., 1983, MNRAS, 202, 995 

Jerjen, E., Bingelli, B. 1997, in the Nature of Elliptical Galaxies; 

The Second Stromlo Symposium, ASP Conf. Ser., 116, 239 
Khosroshahi H.G., Wadadekar Y., Kembhavi A., 2000, ApJ, 533, 

162 

Lima Neto G. B., Gerbal D. & Marquez I., 1999, MNRAS, 309, 
481 

Mazure A. & Capelato H.V. 2002, A&A, in press 

MoUenhoff C, Heidt J., 2001, A&A, 368, 16 

Moriondo G., Giovanardi C, Hunt L.K., 1998, A&AS, 130, 81 

Prieto, M., Aguerri, J.A.L., Varela, A.M. & Munoz-Tunon, C, 

2001, A&A, 367, 405 
Sandage, A., 1961. The Hubble atlas of galaxies, Washington: 

Carnegie Institution 



Seigar, M. & James, P., 1998, MNRAS, 299, 672 
Sersic J. 1968. Atlas de Galaxias Australes Cordoba: Obs. As- 
tronomico 

Simien, F., de Vaucouleurs, G., 1986, ApJ, 302, 564 
Simonneau E., Varela A., Munoz-Tunoz C, 1998, II Nuovo Ci- 
mento, 113B, 927 



Simonneau E., Prada F., 1999, (astro-ph/9906151) 
Stark A., 1977, ApJ, 213, 368 

Trujillo, I., Aguerri, J.A.L., Cepa, J., Gutierrez, CM., 2001, MN- 
RAS, 321, 269 

Trujillo, I., Graham, A.W., Caon, N., 2001, MNRAS, 326, 869 
Varela A., Muiioz-Tuiioz C, Simonneau E., 1996, A&A, 306, 381 
Young, C.K., Currie, M.J., 1994, MNRAS, 268, Lll 
Young, C.K., Currie, M.J., 1995, MNRAS, 273, 1141 
Zhao H., 1997, MNRAS, 287, 525 



APPENDIX A: MASS DENSITY 
APPROXIMATION PARAMETERS 

The table Al shows the values of the parameters that appear 
in the mass density approximation (Eq. 7). 



APPENDIX B: GRAVITATIONAL POTENTIAL 
AND FORCES OF A TRIAXIAL HOMOLOGOUS 
STRUCTURE 

Bl Gravitational Potential 

The gravitational potential at position x = [x, y, z) may be 
written as 



<l>(x) = TvGabc 



lo ^/{r + a^){T + b^){r + c'') 
(Chandrasekhar 1969, p. 52, theorem 12), with 



(Bl) 



(B2) 



It follows from Eq. (Bl) that the potential at an internal 
point is a result of two contributions: that due to the ellip- 
soid interior to the point x considered, and that due to the 
homoeoidal shell exterior to x: 

$(x) = 4TTGaf3 ^ 



F ( arcsin \/T^^^, 
C p(C )F arcsin 



1-/32 



{1-I3^)C^ 
C'2 + A ' 



1-/32 



with F(p,q) the Elliptic integral of the first kind 
(Abramowitz & Stegun 1964, p. 589), and with the restric- 
tion 

„2 



+ 



y 



C'2-hA ' (C'a)2 + A (C'/3)2-f A ^^"^^ 

The strength of the central potential decreases as the 
triaxiality of the object increases: 



$(0)^ = -l-^F { arcsin 



1 - q2 
1-/32 



(B5) 
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Table Al. Parameter values of the mass density approximation 



n 


y 


P 


hi 


h2 


h3 


Max. Rcl. Error (%) 


n ^ 

U.tJ 


-u.ouuuu 


1 nnnnn 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


n nnn 


1 n 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


n nnnnn 
u.uuuuu 


nnn 


1 ^ 




n (\^ nn7 

U.UJ-UU / 


n n79'^7 

-U.U 1 I 


n 9nn4s 

-U.ZUU4:0 


n m fiA7 


1 on 


9 n 


U.^ 1 1 1 o 


n 77zlQ1 
U. ( 1 4:y ± 


n nAQfi'^ 


-U. lOUtJU 


n ns9S4 








n Riin7i 


n n'^'^1 

-u.uooio 


n 1 9n7n 

-U. -LZU ( U 


n 1 A'tnn 
u. i^oyu 


n 070 

U.U 1 u 


o.u 




u.o / uoy 


n n99^i9 

-u. UZZOZ 


n noRi 1 


n 1 Qfisn 

u. ±yuou 


n ofio 
u.uuu 


O.O 


n /loocn 
U.4yzoU 


n Qfim A 

u.oyyi4 


-U.U1d4o 


-U.U ( yiy 


U.z4iuo 


U.UoU 


4.0 


0.50325 


0.91365 


-0.01248 


-0.06747 


0.27969 


0.020 


4.5 


0.51140 


0.92449 


-0.00970 


-0.05829 


0.31280 


0.020 


5.0 


0.52169 


0.93279 


-0.00773 


-0.05106 


0.34181 


0.015 


6.0 


0.55823 


0.94451 


-0.00522 


-0.04060 


0.39002 


0.005 


7.0 


0.58086 


0.95289 


-0.00369 


-0.03311 


0.42942 


0.005 


8.0 


0.60463 


0.95904 


-0.00272 


-0.02768 


0.46208 


0.004 


9.0 


0.61483 


0.96385 


-0.00206 


-0.02353 


0.48997 


0.004 


10.0 


0.66995 


0.96731 


-0.00164 


-0.02053 


0.51325 


0.005 



As reported in Ciotti (1991), the models with low n 
have an inner (r < r^) potential which is much flatter than 
models with high n. As the triaxiality increases there is no 
important change to the gradient of the gravitational po- 
tential along the semimajor-axis; the main effect is to shift 
the gravitational potential profile inwards from the spherical 
case, resulting in a lower potential at intermediate radii. 

B2 Gravitational Forces 

The gravitational forces for a triaxial structure are given by 
the expression: 

axi 

f 1 CV(C) rfC' 

Jo (C'««/a)' + A ^(^'2 + A)((aC')2 + A)((/3C')' + /(^ ^0 ' 
i = 1,2,3 (B6) 
with ai = a, a2 = b, 03 = c, and 

,ii;^,3 [(C a^/a)^ + X] 
The restrictions given in Eqs. (6) and (B4) also apply here. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science I^T^^X style file. 
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